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be  used  to  make  terrain  displays  with  different  characteristics.  Emphasis 
is  on  the  2-photograph  stereoscopic  hologram  which  is  called  a "holographic 
stereomodel".  Techniques  to  enhance  certain  characteristics  of  holographic 
displays  such  as  color  rendition  and  efficient  use  of  illumination  are 
examined  and  possible  uses  of  holography  in  tasks  related  to  map  making 
are  suggested. 
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units.  However,  the  dialable  voltage  cutoff  in  the  battery  charger  will 
permit  recharge  of  silver-cadmium  and  lead-acid  batteries,  among  others. 
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to  the  requirements  of  digital  terrain  data  collection.  Therefore,  this 
approach  is  being  investigated  further  for  implementation  in  an  array 
of  fast,  micro-programmable  processors  to  provide  a benchmark  of  matching 
system  parameters  and  performance. 


14 


UNCLASSIFIED 


security  classification  of  this  page  (When  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  report  number  2.  GOVT  ACCESSION  NO. 

ETL-0091  AD-A035  156 

3.  RECIPIENT'S  CAT  ALOG  NUMBER 

4.  title  (and  SuStirl*; 

DIGITAL  CARTOGRAPHIC  STUDY  AND  BENCHMARK 

SECOND  INTERIM  TECHNICAL  REPORT 

5.  Type  of  report  a period  covered 

Contract 

6.  performing  org.  report  number 

7.  AUTHORf*; 

D.J.  Panton 

M.E.  Murphy 

a.  contract  or  grant  NUMBER(a> 

DAAG53-75-C-0195 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Control  Data  Corporation 

2800  East  Old  Shakopee  Road 

Minnpapnlis,  Minnp<;nta  88440 

10.  program  element,  project,  task 

AREA  A WORK  UNIT  NUMBERS 

n.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.S.  Army  Engineer  Topographic  Laboratories 

Fort  Belvoir,  Virginia 

12.  REPORT  DATE  | 

December  1975 

13.  NUMBER  OF  PAGES 

28 

U.  MONITORING  AGENCY  NAME  A ADORESS('f(  diUereni  from  Controlling  0/(lce; 

15.  SECURITY  CLASS,  (o/ th/«  report; 

Unclassified 

ISA.  OECLASSf  FI  CATION/ DOWN  GRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thim  Report) 

Approved  for  Public  Release;  Distribution  Unlimited 

17.  DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  In  Block  20,  H different  from  Report) 

IS.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Condnua  on  rereree  aide  It  necememry  and  Idantlfy  by  block  number) 

Stereo  Image  Matching  Microprocessors 

Altorithm  Distribution  Match  Point  Reliability 

Block  Matching  Algorithm  Correlation  Patch  Shaping 

Microprogrammable  Processors  Flexible  Processor{FP)  Array 

20.  ABSTRACT  (Caatkaue  marereten  ai^  M neeeaeaey  aad  fdietUlty  by  block  number) 

This  report  describes  the  practicality  of  implementing  the  automatic  process 
of  stereo  image  matching  on  a configuration  of  extremely  fast  micro- 
programmable  processors.  The  report  includes  a three-phase  program  designed 
to  bridge  the  gap  between  mathematical  algorithm  development  and  an  actual 
hardware  benchmark  realization.  The  report  concludes  with  the  block  matchini 
algorithm  being  logically  reconstructed  and  ready  for  microprogramming. 

DO  /.STt,  I4n  eiMTIOH  OF  I HOV  S>  IS  OBSOLETE  15  UNCLASSIFIED 


SECUBITY  CLSiStFICATlOH  OF  THIS  PAGE  (•■hwt  Dmim  EnlarMf) 

— - 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  fWhtn  Pmlt  F.nffd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

Y.  REPORT  NUMBER  | 

ETL-0092  i 

2.  GOVT  ACCESSION  NO. 

AD-A035  157 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (trd  SubtItU) 

5 TYPE  OF  REPORT  & PERIOD  COVERED 

DIGITAL  CARTOGRAPHIC  STUDY  AND  BENCHMARK 

Contract 

3RD  INTERIM  TECHNICAL  REPORT 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^*; 

8.  CONTRACT  OR  GRANT  NUMBERr«J 

D.  J.  Panton 

DAAG53-75-C-0195 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Control  Data  Corporation 

2800  East  Old  Shakopee  Road 
Minneanol is.  Minnesota SS44Q 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  REPORT  DATE 

September  1976 

U.  S.  Army  Engineer  Topographic 

Laboratories 

IS.  NUMBER  OF  PAGES 

38 

14.  MONITORING  AGENCY  NAM^ft  ADDRESS^^ff  (f///•r«nr  from  Controlling  Oitico) 

IS.  SECURITY  CLASS,  fpf  Ififa  r«pe>ft; 

Unclassif ied 

1S«.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thl»  Roport) 

Approved  For  Public  Release;  Distribution  Unlimited 

17.  DISTRIBUTION  STATEMENT  (of  tho  obotrmet  mtfrod  In  Block  20,  If  difforont  from  Roport) 

16.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continuo  on  rororoo  oldm  If  uoeomoory  mnd  Idontity  by  block  number) 

Benchmark  algorithm 

Correlation  patch 

Parallel  Arrays 

Ml  abstract  fCmiflMiio  em  roeotwo  oM*  H nmommmoey  mW  Identity  by  block  rwmbee) 

This  report  describes  the  work  and  findings  of  Phase  C,  which  Is  to  Implement 
and  test  on  fast,  mlcroprograiitmable  processors,  the  stereo  matching  algorithm 
that  was  designed  and  analyzed  under  Phases  A and  B.  This  Implementation  was 
performed  in  terms  of  a benchmark  to  shed  some  light  on  the  practicality  of 
such  a concept  and  to  uncover  the  advantages  and  disadvantages  of  this 
Particular  kind  of  parallel  processing  application.  The  results  stem  from 
actually  performing  the  benchmark  rather  than  from  a paper  study  outlining  Cie 

W)  » jSSw  M73  Eoi-noti  of  » MOV  6»  is  cmmouete  unclassified 


SECUMTV  CLASSIFICATION  OF  THIS  PACE  (Whmt  DM*  Bnitnd) 


20. 


continued 


and  accuracy.  It  was  found  that  it  takes  about  11  hours  of  GDC  6600 
processing  to  complete  the  correlation  and  matching  of  a typical  9X9 
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lay  the  groundwork  for  handling  non-central  perspective  photography. 
Previous  developments  and  algorithm  logic  modifications  have  been  re- 
ported in  a rather  piecemeal  fashion  over  the  first  three  Phases.  This 
report  combines  all  these  developments  into  a consistent  description  of 
the  matching  algorithm  as  it  appears  to  date,  including  the  modifications 
of  Phase  D. 
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expect  tolerable  y-parallaxes.  Therefore,  either  the  transparencies  must 
be  tilted,  or  the  HOC  should  operate  essentially  as  a comparator,  with 
correlation  taking  place  in  the  image  space.  For  the  latter  case,  one 
transparency  is  fixed  and  the  other  is  moved  in  steps  in  both  the  x-  and  y- 
di recti ons.  The  result  is  that  for  each  array  element  one  gets  x,  y 
coordinates  on  one  transparency  and  p^,  Py  parallaxes  from  the  other. 

The  results  of  two  preliminary  tests  using  aerial  photography  are  given. 
They  indicate  accuracies  of  about  one  to  two  parallax  steps  (one  step  = 

25  urn  in  the  tests).  However,  the  tests  are  so  limited  that  the  results 
should  be  considered  only  as  preliminary,  indicating  the  feasibility  of 
the  system  as  an  area  correlator.  Further  development  and  more  extensive 
testing  are  recommended  before  the  potential  and  accuracy  of  the  HOC  system 
may  be  detennined. 
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map  and  chart  requirements:  0 - Not  detectable,  the  map  element  cannot  be 
discerned  or  located  from  either  type  of  LANDSAT  data;  1 - Detectable,  map 
element  can  be  detected  but  not  identified  from  the  type  of  LANDSAT  data 
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indicated;  2 - Identifiable,  map  element  can  be  detected  and  recognized 
as  a particular  type  of  feature  from  the  LANDSAT  data  indicated,  e.g. 
road,  canal,  etc.,  collateral  information  may  be  required  to  reach  this 
analysis  level;  3 - Classifiable,  LANDSAT  data,  with  the  use  of  all 
available  collateral  information,  can  provide  the  information  required 
for  the  map  element  including  all  required  measurements,  e.g.  width,  length, 
and  areas.  It  was  concluded  that  LANDSAT  1,  2,  3 MSS  data  is  compatible 
with  National  Map  Accuracy  Standards  and  can  be  used  to  update  the  map 
elements  on  map  scales  1:1,000,000  through  1:250,000,  although  many  of 
the  cultural,  hydrographic,  and  botanical  elements  may  be  unclassifiable. 

The  improved  systems  capabilities  of  LANDSAT  4 may  provide  a method  for 
updating  map  scales  1:1,000,000  through  1:50,000.  However,  many  of  the 
required  cultural  and  hydrographic  map  elements  may  remain  unclassifiable 
even  with  the  Thematic  Mapper  system. 
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Insight  is  given  into  the  nature  of  depolarization,  but  explicit  results 
for  the  depolarized  terms  are  not  obtained  at  this  time  because  of  the 
complexity  and  difficulty  of  the  solution.  Some  of  the  conclusions  of 
this  work  are 

1.  k theory  has  been  developed  for  computing  the  like  polarized 

(HH  and  VV)  radar  backscatter  coefficients  from  certain  types  of  vegeta- 
tion by  using  a vector  renormalization  approach. 

2.  No  rigorous  quantitative  comparison  of  theory  with  experiment 
was  possible;  however,  qualitative  comparisons  indicate  reasonable  agree- 
ment . 


3.  Although  no  explicit  solution  was  obtained  for  the  depolarization 
components,  it  was  learned  that  one  cause  of  depolarization  is  the  aniso- 
tropy associated  with  the  correlation  function  of  the  dielectric  fluct- 
uations . 
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coordinates  into  their  correct  ground  positions  with  relief  displacement 
errors  removed.  The  second  system  uses  the  mono  solution  for  mono  input 
imagery  and  a full  analytical  mapping  solution  for  stereo  input  imagery. 
The  systems'  design  evolution,  theory  of  operations,  technical  character- 
istics, estimated  costs,  field  deployed  van  layout,  and  cost  effectiveness 
comparisons  are  presented. 
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the  principle  of  array  algebra  in  a piecewi.se  translocation  algorithm 
is  applied.  In  this  approach  the  non-gridded  epipolar  coordinates  are  first 
converted  to  regularly  spaced  elevation  data  and  then  subsequently  com- 
pacted using  the  methods  of  array  prediction.  In  addition  to  analyzing 
the  mathematical  equations  required  for  terrain  compaction,  the  com- 
putational requirements  were  analyzed  for  both  sequential  and  parallel 
processors . 
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RGSS  performance  in  open  traverse  is  unlikely  without  vertical  deflection 
end  calibration.  In  addition,  the  keynote  of  successful  RGSS /Grad iome ter 
integration  will  be  control  and  compensation  of  system  bias  and  low 
frequency  error  sources. 
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20.  continued 


erometer  Errors. 

This  report  concluded  that  (1)  the  gyro  and  accelerometer  errors  limit 
ability  to  derive  deflections  from  frequent  zero  velocity  updates  (at 
vehicle  stops);  gradiometer  vs  improved  inertial  system  components  are 
trade-off  possibilities,  (2)  post-mission  deflection  recovery  to  0.5 
sec  appears  feasible  with  current  IPS  hardware  augmented  by  high  quality 
deflection  values  at  both  traverse  end  points  and  a moving-base  gradio- 
meter with  field  performance  comparable  to  that  recently  demonstrated  by 
laboratory  prototypes  (RMS  noise  < 10  EU,  10  second  average,  (3)  Real- 
time deflection  recovery  to  0.5  sec  is  ultimately  possible  in  an  en- 
vironment of  moving-base  gradiometer  availability,  highly  accurate  end 
point  deflection  data  and  moderately  improved  RGSS  inertial  components, 
(4)  the  dominant  factor  limiting  RGSS  deflection  recovery  both  with  and 
without  gradiometer  aiding  is  the  quality  of  vettical  deflection  data 
which  is  independently  available  at  traverse  end  points. 
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sate  1 1 ite- to-sate 1 1 ite . Each  Geoceiver  type  may  use  either  a fixed  or  a 
variable  averaging  time.  Included  in  the  implementation  of  the  Geoceiver 
measurements  is  the  Hopfield  model  for  computing  trpospheric  refraction 
correct  ions . 

The  above  modifications  have  been  implemented  both  on  the  L^'IVAC 
1108  and  the  CDC  6400  versions  of  Fotonap. 
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arrays.  Since  each  array  produces  four  simultaneous  video  signals, 
the  system  could  be  utilized  as  an  input  device  for  a parallel  processor 
such  as  the  STARAN.  The  video  signals  can  be  multiplexed  onto  a single 
line  if  desired. 
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Electron  Beam  Recording  (EBR)  was  studied  in  the  context  of  its  potential 
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several  such  materials  were  found  to  be  suitable. 
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A program  of  periodic  inspection  and  continuing  evaluation  is  required  to 
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mission  --  as  well  as  the  data  on  the  IPS  velocity  errors  at  each  stop. 

The  purpose  of  this  report  is  to  develop  the  equations  for  the  position, 
velocity,  and  tilt  angle  errors  into  a useable  algorithm  for  the  optimal 
estimation  of  the  deflections  of  the  vertical.  As  a result,  three 
major  analytic  tasks  are  presented.  This  mathematical  development  enables 
the  production  of  a machine  algorithm  for  use  with  actual  data.  The  final 
section  of  the  report  contains  a coded  version  of  the  algorithm  with 
explanatory  comments. 
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year  field  program  sponsored  by  the  U.S.  Army  Natick  Laboratories  and 
included  ground  station  temperature,  relative  humidity,  precipitation, 
evaporation,  hours  of  sunshine,  and  tower  measurements  of  temperature,  dew 
point  temperature,  and  wind  speed,  incoming  and  reflected  solar  radiation, 
incoming  and  outgoing  infra-re^  radiation  and  sub-surface  temperature 
profiles  in  the  forest  and  in  a cleared  area  within  the  forest  region. 

The  results  discussed  here  include  a climatological  survey  of  the  e.x- 
perimental  iste  as  well  as  a comparison  of  the  monsoonal  cycles  at  a 
number  of  widely  dispersed  sites  in  Southeast  Asia;  an  analysis  of  the 
piofiles  of  temperature  and  wind  speed  in  the  forest  and  cleared  area; 
diurnal  variations  of  static  stability  and  Richardson  number  repre- 
sentative of  several  different  periods  of  the  year;  a summary  of 
calculations  of  the  roughness  characteristics  of  the  forest  canopy  and 
the  variability  thereof;  spectral  analyses  of  kinetic  energy  above  and 
within  the  forest  canopy  during  both  monsoonal  flow  regimes;  estimates 
of  the  albedo  of  the  for  st  and  the  net  fluxes  of  radiational  energy  for 
several  periods  of  the  year;  estimates  of  sensible  and  latent  heat 
fluxes  above  the  forest  canopy. 
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